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The A + B diffusion-limited reaction, with equal initial densities a(0) = 6(0) = no, is studied by 
means of a field-theoretic renormalization group formulation of the problem. For dimension d > 2 
an effective theory is derived, from which the density and correlation functions can be calculated. 
We find the density decays in time as a, 6 ~ C^fK{Dt)~'''^'^ for d < 4, with A = no — C'n^^^ + ■ ■ 
where C is a universal constant, and C' is non-universal. The calculation is extended to the case 
of unequal diffusion constants Da Db, resulting in a new amplitude but the same exponent. For 
d < 2 a controlled calculation is not possible, but a heuristic argument is presented that the results 
above give at least the leading term in an e = 2 — d expansion. Finally, we address reaction zones 
formed in the steady-state by opposing currents of A and B particles, and derive scaling properties. 
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I. INTRODUCTION 

Diffusion-limited chemical reactions are known in lower 
dimensions to exhibit anomalous kinetics [Q,^- That is, 
the evolution of the density depends strongly on fluctua- 
tions, and cannot be derived from mean-field rate equa- 
tions. In this paper we apply renormalization group (RG) 
techniques to the two-species reaction A + B ^ %, with 
the goal of determining systematically the effects of these 
fluctuations. 

The model for the A + B ^ % reaction involves 
two types of particles, both undergoing diffusive random 
walks, and reacting upon contact to form an inert parti- 
cle. In the density rate equation approach it is assumed 
that the A and B particles densities a and b are uniform, 
and that reactions occur at a rate proportional to the 
product ab, giving 
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with rate constant F. In the case of equal initial densities 
a(0) = 6(0) = no, the solution goes as a, 6 ~ {Tt)~^ 
asymptotically, with an amplitude which is independent 
of the initial density. 

It was first suggested by Ovchinnikov and Zeldovich 
1^ , and later demonstrated by Toussaint and Wilczek , 
that relaxing the assumption of uniformity yields a slower 
density decay. In particular, Toussaint and Wilczek made 
the observation that if the two species have the same 
diffusion constants Da = Db = D, then the density 
difference a — b obeys the diffusion equation. As a result 
they found, by using central limit arguments to calculate 
the fluctuations in a — 6 due to equal density, random 



initial conditions, the asymptotic densityB 
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where d is the dimension of space. Comparing with the 
result of the rate equation, we see that for d < 4 the 
asymptotically dominant process is the diffusive decay of 
the fluctuations in the initial conditions. 

Using a particular version of the model, Bramson and 
Lebowitz confirmed rigorously the decay exponent of 
Toussaint and Wilczek, finding for d < 4 



(1.3) 



where Cd is some constant which depends on the dimen- 
sion d 1^,^. In their treatment they demonstrated that 
the two species are asymptotically segregated for d < 4. 
This seg regation was assumed in reference Q in deriving 
Eq. (|l|). 

Numerical simulations have confirmed the value — d/4 
for the decay exponent in one QJ^, two 1^,^, and three 
dimensions For all of these simulations restric- 

tions were placed on the occupation number per site, 
and usually only single occupancy allowed. In the one- 
dimensional simulation of Toussaint and Wilczek the ini- 
tial density was varied, and reason abl e agreement was 
found with their analytic result, Eq. (1.2) However, in 
higher dimensions the ^/uq amplitude dependence, when 
tested, has not been observed [^,1). In the former case 



* There is a misprint in the d — 3 amplitude of Ref. 
(19c). 
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the initial average occupation number per site was kept 
low, whereas for the higher dimensional simulations it 
was necessary to start with a nearly full lattice in order 
t o re ach the asymptotic regime. This suggests that Eq. 
(1.2) might not be a universal result, but rather a limit 
for small initial density hq. 

While (i = 4 appears to be the upper critical dimen- 
sion for homogeneous initial conditions, this is not the 
case when the two species are initially segregated, where 
instead the upper critical dimension is found to be c? = 2 
pO| . That is, as a result of the segregation, a localized 
region forms in which nearly all reactions occur. This re- 
action zone exhibits scaling behavior, and the character- 
istic exponents are independent of the dimension d when 
d > 2, but crossover to dimension-dependent values for 
d < 2. Hence, one of our goals in applying RG techniques 
is to better understand the role of the dimensions d — 2 
and d = 4. 

The problem can be mapped to a field theory by start- 
ing from a master equation description of the model [ pT| . 
From an analysis of the field theory we find that there is 
an upper critical dimension d^ — 2, which is associated 
with the stochastic processes of reaction and diffusion. 
Hence, for d > 2 one can replace the full field theory 
with an effective theory, which is valid for asymptotically 
late times, while for d < 2 one must instead perform an 
explicit renormalization group calculation. 

The effective theory for d > 2 is equivalent to the de- 
terministic partial differential equations 

5ta = DaV^o - Tab dtb ^ DsV^b - Tab (1.4) 



with stochastic, non-negative effective initial conditions. 
In deriving the effective theory we find that the initial 
distribution is finitely renormalized due to the presence 
of relevant initial terms, the analog of surface terms for 
a t = boundary in a d -I- 1 dimensional theory. The re- 
sulting distribution can be characterized by a parameter 
A which depends nonuniversally on the initial density. 

We demonstrate explicitly that from these equations 
follows generally the asymptotic segregation of the A, 
B particles when d < 4, and subsequently the universal 
decay exponent — d/4. However, the amplitude of the 
density decay depends on the initial conditions, and is 
therefore nonuniversal. It is important to note that if 
one uses instead central limit arguments to calculate the 
initial distribution which should be fed into ( |1.4D , then 
one is implicitly making the assumption that these equa- 
tions hold for all times, rather than just asymptotically. 
Such an assumption will get the exponent correct, but 
we claim that it does not, in general, predict the cor- 
rect amplitude because it neglects the dynamics at short 
times. 

For 2 < d < 4 and Da = Db wc find 
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where the angular brackets denote averages over both the 



processes of reaction and diffusion and the initial condi- 
tions. Here A is the coupling constant of the induced 
initial terms, and can be calculated as an expansion in 
the initial density, giving 
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where AoS is a nonu niversal effect ive ra te constant, de- 
fined in section II C and used in ( 2.2S ). Hence, in the 
small no limit the amplitude is universal, and we re- 
cover the result of Toussaint and Wilczek, Eq. ( |l.2| ). The 
higher order terms in no are nonuniversal, and offer a 
possible explanation for the deviation from y^no behav- 
ior found in the siin ulatio ns |p| ,^ . 

Our results ( |l.5| ) and (1.6) appear to disagree with 
those of Bramson and Lebowitz Q , for d = 4 as well as in 
the case above. However, we stress that since the depen- 
dence of the amplitude on the initial density is nonuni- 
versal there is no explicit contradiction. Our model is 
defined by a continuous time master equation in which 
the reaction occurs at a rate A, and multiple occupancy 
per site of each particle type is allowed. Bramson and 
Lebowitz also study a continuous time model with mul- 
tiple occupancy allowed, but with an instantaneous re- 
action ||,^. In this case a lattice site can only contain 
one type of particle. We use a finite reaction rate since 
this is convenient for mapping to the field theory, and 
because it allows one to determine better the extent of 
universality. However, we cannot directly relate our re- 
sults to those of Bramson and Lebowitz, since the field 
theory techniques we use are no longer valid in the limit 
A oo, to which their model corresponds. We note that 
if our results should be valid for large but finite A, then 
Acff, given by ( |l.6D , goes to a limiting value of the order 
/i''~2, where h is the short distance cutoff. 

For d < 2 the full field theory and the subsequent 
renormalization must be considered. We find that the 
field theory may be exactly renormalized, as was shown 
previously by Peliti for the one-species reactions A+A^ 
A and A + A ^ |jl^. However, the e-expansion calcula- 
tion of observables requires non-perturbative sums over 
all orders of the initial density no and the parameter A, 
and while these may be carried out straightforwardly in 
the one-species reaction [|3| , we are unable to apply these 
methods to the present case beyond the leading order in 
e = 2 — d. Thus, at least in this approach, we are not yet 
able to establish even the power law for the density 
to all orders in e, although we believe it to be true. 

We also consider the case of unequal diffusion con- 
stants. Da 7^ Db, and show from the effective theory 
that in the small no limit 



{a{S)) ^ VQR^ («(^ = 0)) 
where S = {Da - Db)/{Da + Db), and 
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Therefore this falls into the same universality class, with 
regard to the decay exponent, as the symmetric case. 

From the effective field theory for 2 < d < 4 it follows 
that the density difference a — & is at late times a gaus- 
sian random field. This, combined with the asymptotic 
segregation a + b ^ W — b\ allows one to calculate any 
correlation function. We calculate exactly the equal time 
two-point correlation functions (a(r)a(O)) and (a(r)6(0)). 

The final topic we discuss is that of reaction zones, 
which form whenever A and B particles are segregated. 
One example of a reaction zone is that which results 
from opposing currents of A and B particles. We ap- 
ply RG methods to this steady-state case, and show that 
the densities and the rate of reaction have universal scal- 
ing forms. The upper critical dimension for this system 
is dc = 2. These results can be extended to apply to 
reaction zones in initially segregated systems pO| , p^ , p^ , 
and also homogeneous systems for d < 4 . 



II. THE MODEL AND THE CORRESPONDING 
FIELD THEORY 

The model is defined by a continuous time master 
equation for the probabiHty P({m}, {n}, t). The set {m} 
denotes the occupation numbers of A particles on each 
lattice site, {n} the occupation numbers of B particles, 
and P is the probability of a given configuration occur- 
ring at time t. The master equation for P reads 

d 
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(2.1) 

where Da, Db are the diffusion constants for A and B 
particles, h is the size of the hypercubic lattice, and A is 
the microscopic reaction rate constant. In the first two 
curly bracket terms, which describe the diffusion of A and 
B particles, i is summed over all sites, and j is summed 
over nearest neighbors to i. 

The initial conditions for P are given by a Poissonian 
distribution, with the average occupation number per lat- 
tice site equal to no for each species. That is. 
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A. Mapping to Field Theory 

The first step in mapping the master equation to a 
field theory is to recast it in a 'second quantized' form, 
following a procedure developed by Doi flj . Two sets of 
creation and annihilation operators — a, a 'for A particles 
and 6, 6^ for B particles — are introduced at each lattice 
site. These obey the usual commutation relations: 



(2.3) 



with all other commutators zero. The vacuum ket is de- 
fined by ai|0) =0 and bi\0) = for all i. In terms of these 
operators the state of the system at time t is defined to 
be 

E Pi{rn},{n},t)l[{air{^r\0). i2A) 

{m},{n} i 

The master equation can be rewritten in terms of this 
state as 



d_ 

di 



\m = H\m, 



(2.5) 



with the operator 

<ij> ^ 

+AE(aI^l-l)a.&.- (2.6) 



(2.7) 



The formal solution of Eq. ( |2.5| ) is 

\m)^e-f'\m)- 



The density and other averages, which are defined in 
the original occupation number representation, can be 
calculated from \4i{t)) by introducing the projection state 



(2.8) 



in terms of which the average is 

mt)))^ E A({m},{n})p({m],{n}, 
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The operator analog A can be derived for any 
A({m}, {n}) by Taylor expanding the latter with respect 
to mi,ni, and then substituting 
Note that 



aldi, fii -> bl'bi 



= ( 16 = 



(2.10) 



for all i, implying that A can be expressed solely in terms 
of annihilation operators by first writing it in normal or- 
dered form. 
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The second quantized version of the model is mapped 
to a field theory by the use of the coherent state rep- 
resentation 18|. The time evolution operator in Eq. 
(2.7) is rewritten via the Trotter formula 



exp(-i/t) = lim (1 - HAt) 

At^O 



t/At 



(2.11) 



The right-hand side, before the At ^ limit is taken, 
can be regarded as a factorization of the operator into 
time slices, and a complete set states inserted between 
each factor. The identity is given in the coherent state 
basis by 



(2.12) 



where \z) is the normalized eigenstate of the annihilation 
operator with complex eigenvalue z: 



|0). 



(2.13) 



Equation (2.12) is generalized to a product over all lat- 



tice sites and particle species, and then the operator 
(1 — HAt) is evaluated between successive time slices, 
r esu lting as A< ^ in a path integral representation of 

(H. 

The corresponding action is 



S 



dtl a 



Da 
D 



V 



+b*(^dt-^S/^y-Xoil-a*b*)ab^ 



-?ioa*(0) - nob*{0) ~ a{tf) ~ b{tf) 



(2.14) 



where the fields a,a*,b,b* originate from the complex 
variables z for each particle type. Time has been rescaled 
by the average diffusion constant D = {Da + Db)/2, and 
the coupling constant is given by Aq = Xh'^/D. The time 
derivatives above come from the overlap between time 
slices, and the other the curly brace terms result from 
the operator H . The remaining terms are not integrated 
over time and represent the random initi al st ate, with 
no = fio/h'^, and the final projection state (2. 



Averages, as defined in ( |2.9| ), are given in terms of this 
action by 



mt))) = 
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where the script V denotes functional integration. The 
functional A[a, b] is found by directly substituting the 
fields a, b for the annihilation operators d, b in A. 

The time i/ of the projection state is arbitrary as long 
tf > t, where t is the time argument of the observable. 
This follows directly from the condition = for prob- 
ability conservation. The final terms can be eliminated 



by making the field shifts a* 
Then the reaction terms are 



1 + a and b* 



-Ao(l -a*6*)a6 



Xo{d + b)ab + Xodbab. (2.16) 



Since the conserved mode a — b plays an i mport ant role 
in the dynamics, it is useful to transform ( 2.14 ) to the 
fields (f>, 4>, ip, 'ip defined by 

a + b - a + b 



V2 



V2 



a — b 



V2 • 
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The_\/2 factors are included so that the derivative terms 
in (2.14) maintain a coefficient of unity. The subsequent 
action is 



5 = j d'^x 

-50VV- Ai 



(2.18) 



where 5 — {Da — Db)/{Da + Db), the couplings are 
•^1 = Xq/^/2 and A2 = Ao/4, and the initial density is 
= v^no- This action is the starting point for our 
analysis. Note that, since we are considering only equal 
initial densities, (a) = (6) = {(f)) /y/2. 

The mapping outlined above is a general technique, 
which may be applied to many different reactions, for 
example, the general two-species annihilation reaction 
mA + nB — > 0. However, in the present work we restrict 
ourselves to the case m = n— 1. 



B. Diagrams and Power Counting 

A perturbatio n expa nsion for a given observable can be 
developed from ( 2.15 ) and ( ^.18 ), and expressed in the 
usual diagrammatic fashion. Th e prop agators for (f>, ^ are 
given by the first two terms in (2.18), and are the diffu- 
sion equation Green's function: G^^{k,t) = G^^{k,t) = 

6"'=^* when t> and G^x = G,i„t, = for t < 0. These 



propagators are represented by solid and dashed lines re- 
spectively. The three- and four-point vertices, which cor- 
respond to the annihilation reaction, are shown in Fig. 0. 
When 6^0 there are also two-point vertices which con- 
nect a (f) propagator to a ip propagator, and vice versa. 
These vertices are wave number dependent, with magni- 
tude fc^. 

In addition there is a source term g^<t><l>it=o) ^ gy Baylor 
expanding the exponential an expansion in powers of 
is generated, where the diagrams giving the coefficient 
have a source of i (j) propagators at t = 0. It is useful to 
introduce the classical density and the classical response 
function, which both involve sums over all powers of n^. 
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FIG. 3. The classical response function, shown as a heavy 
line, is given by the sum of all possible tree diagrams con- 
nected to a single propagator. 



FIG. 1. Pro pagat ors and vertices for the full theory, given 
by the action (2.18). 



(a) 




(b) 



FIG. 2. The classical density, represented by a wavy line, 
is given by (a) the sum over tree diagrams and (b) an integral 
equation. 



These quantities are important because it is found that 
under renormahzation flows to a strong coupling limit, 
and these sums are still meaningful in this limit |13|| . The 
term 'classical' refers to absence of loops in the diagrams. 

The classical density is defined to be the sum of all 
tree diagrams which contribute to the average (0), as 
shown in Fig. ||. Note that these diagrams contain only 
propagators, because of the three-point vertices in (2.18). 
This sum obeys an integral equation which can be solved 
exactly, giving 



1 + n^Xit 



(2.19) 



The classical response function is defined to be the (/) 
propagator with all possible tree diagrams branching off 
to t = 0, as shown in Fig. Again this can be solved 
exactly, giving 



(0(k, i2)0(-k,ti))ci = e 



1 -I- n^\it2 
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In order to renormalize the field theory we must first 
determine the primitive divergences, for which we con- 
sider the following dimensional analysis. There is a rigid 



constraint that [^0] = [^tp] = ^ , where k has the dimen- 
sions of wave number. If we take the dimensions of the 
conjugate fields to be [01 = [t/S] = fc", as was done for the 
one-species reaction |1^, then a general vertex ■ip'^ (jp ip^ ip^ 
is found to be relevant only for k + l < 2 and d < 2. 
Next, we observe that it is not possible to genera te an y 
vertices with fc = 1 or £ = 1 from the vertices in (2.18). 
Th erefor e all relevant vertices are exactly those present 
in ( 2.18D , with an upper critical dimension dc — 2. We 
will discuss the renormahzation of the theory in section 



[lie, and for now focus on the case for d > 2. 



To elucidate the crossover which occurs for d > 2 it 
is useful to consider rescaling the fields by dimensionful 
parameters, which is consistent as long as the conjugate 
fields are rescaled accordingly. Under such a rescaling 
the couplings Ai and A2 behave differently, although orig- 
inally they are both proportional to Aq. In particular, we 
can take 



^/\, 



Xi(t> t/j^ip/Xi V 



AiV-, 
(2.21) 



which has the result of setting the 0(0^ — ■0^) coupl ing to 
unity, while leaving A2 unchanged. The rescaling (2.21) 
also results in — > Xin^. This is the proper quantity 
to study when addressing issues of relevance and irrele- 
vance, which can be seen by studying the diagrams gen- 
erated by the action (2. IS): whenever an additional t — 
line is added with weight n^, there is an additional Ai re- 
quired to connect it. In this system of units, then, one 
finds that 



[A2 



r 2-d 



[Xifid 



(2.22) 



Therefore there exists a critical dimension dc — 2, above 
which A2 flows to zero. Doing the complete power count- 
ing method with the rescaled fields yields the same result. 
The initial density is a strongly relevant parameter for all 
d. The diffusion constant difference S is always marginal 
whenever it is not zero. 

Before turning to the consequences of the power count- 
ing for the case d > 2, we mention another approach 
to this problem, which is to integrate out the conjugate 
fields ip and cj) in ( |2.1S| ). This leads to the equations of 
motion (for (5 = 0) 

^0 = V20-Ai02 + AiV/+?70 (2.23) 
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where 77^, ry^ are multiplicative, complex noise terms p9[ . 
It is important to note that the physical density is not the 
field (f), but rather the average of 4> over the noise terms. 
These equations, without the noise terms included, are 
often taken as the starting point for analysis, but this 
approach is not generally valid. As we will show in the 
next section, one can neglect the noise terms only for 
d > 2 and for asymptotically large times. 

Equation (2.24) can be simplified in any dimension, 



since it is a linear equation, by averaging over the noise. 
This is an average over the stochastic process of diffusion, 
and not over the initial conditions. Then the averaged 
field ("0) obeys the simple diffusion equation for any given 
initial configuration. 



C. Effective Field Theory for d > 2 

From the dimensional analysis and power counting 
above it follows that for d > 2 the full theory given by 
(2.18) can be replaced by an effective theory in which 
A2 = and Ai Aoff (Ai, A2, A), where A is a wave 
number cutoff, of the order of the inverse lattice spac- 
ing. However, in constructing such an effective theory 
one has to consider all possible relevant terms, consistent 
with the symmetry of the theory, which might be gener- 
ated under renormalization. In order to identify these 
terms we note that this problem is analogous to that of a 
semi-infinite system in equilibrium statistical mechanics 
in d + 1 dimensions, the analog of the boundary being 
the hyperplane t = 0. While one finds, in the semi- 
infinite equilibrium case, that the bulk critical properties 
do not depend on the surface terms, nonetheless one ex- 
pects surface terms to contribute to correlation functions 
which involve fields on the boundary |Q . All observables 
in our problem are given by such correlation functions, 
since all diagrams originate with the 71^0(0) term. There- 
fore we must check for all relevant initial terms, the t = 
analog of the surface terms, which might be generated, 
as well as those of the bulk. As mentioned above, the 
only relevant bulk term is that of Ai . 

The proper framework for determining which terms 
are relevant is via the rescaled fields ( ^.21 ). Therefore, 
for an init ial t erm of the type (A^™'"-' /m!n!)(^™-0"|i=o 
added to (2.18) we consider the dimensions of the cou- 
pling [A™+'^'"A^'«)] = fc(«+m)(2-d)+d^ rpj^jg pQ^^^gj. 

also follows from calculating the number of vertices re- 
quired to attach a, t = vertex of A^"^'"-'> to a given 
diagram. These terms are relevant when 



d < 



2{n -j- m) 
n + m — 1 



(2.25) 



If m -f n = 1 then the initial term is relevant for all d. 
The case m = I corresponds to the initial density, which 
has already been demonstrated to be relevant. For the 
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FIG. 4. The initial term A is generated by diagrams of 
the form (a). The tree diagrams in (b) give the leading order 
contribution for small no . The leading order corrections come 
from the diagrams (c). For A'-^'"' the same diagrams would 
be used, but with the opposite sign for the A2 vertex on the 
left. 



case n = 1 we first address a symmetry of the theory. 
When starting with equal initial densities the system 
is invariant under exchanging A ^ B and Da Db- 
Therefore the action must be invariant under the trans- 
formation ((/), 0, -!/;, t/i, 5) — > (</), (/), —■(/', —ijj, For what 
follows we will consider only the case (5 = 0, or Da — Db, 
in which case the symmetry forbids the generation of a 
initial term A^"'-*^^-)/). In section IV the case (5^0 will 



be discussed, and it will be demonstrated that again no 
n = 1 initial term is generated. 

For m + n = 2 symmetry allows only the genera- 
tion of A(2^o) and A^O-^). Below we will address the 
calculation of these quantities, and demonstrate that 
/\(o.2) _ _^(2,o) ^ These terms are relevant when- 



ever c? < 4, as can be seen by equation (2.25), and there- 
fore must be considered when constructing an effective 
theory for 2 < d < 4. In fact, it will be shown that the 
term (A/2)'0^ is solely responsible for determining the 
asymptotic decay of the density. This is an important 
point. This system is dominated by initial terms, as op- 
posed to the one-species reaction. Therefore techniques 
which utilize homogeneous source terms and look for a 
bulk steady state will not work for this problem. Since 
this initial term dominates the asymptotic behavior of 
the density, we identify d* = 4 as a second critical di- 
mension of the system. 

Higher order initial terms will also be relevant in the 
range 2 < d < 3. In fact, as d ^ 2 one finds that all 
initial terms become relevant. While this seems to be 
an extreme complication, it is in fact possible to calcu- 
late exactly the asymptotic density for 2 < d < 4 and 
demonstrate that it is independent of such terms. This 
will be presented in the next section. We now turn to 
the calculation of the parameter A. 

The diagrams which must be considered in calculating 
an effective initial term (A/2)-!/'^ are all those in which 
two -i/i lines exit from the left, as shown in Fig. ||(a). 
The sum of these diagrams gives rise to an effective term 
/(t)'0(i)^ in the action. If the function f{t) goes to zero 
for large t, and is sharply peaked enough that Jg°° dtf{t) 
is finite, then a coarse-graining in time gives f{t)'4){ty' 
(A/2)(5(i)?/;(0)^, where both quantities are understood to 
be integrated over t, and A = 2 dt'f{t'). To calculate 
this parameter A we consider first the subset of diagrams 
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given by the tree diagrams, as shown in Fig. ^(b). These 
diagrams sum to give fo{t) — — A2/(l + n^Xit)'^, and so 



An = -2A2 



dt 



1 



2n^ 



(2.26) 



Therefore we conclude that this set of diagrams generates 
an effective initial term Aq — 2n0A2/Ai, or, in terms of 
the parameters in the master equation Aq = no, 
the initial density of each species. This will be shown 
to be the leading order term for a small no expansion of 
A. The width of the function /o(t) is given by (n^Ai)"^, 
and therefore we expect this coarse-grained picture to be 
valid for times t ^ (n^Ai)~^. 

We can group all the diagrams in the full theory ( |2.1S| ) 
which are of the form specified in Fig. ^(a) in the fol- 
lowing way. There is a vertex A2 which is the leftmost 
vertex in the diagram. The lines coming into this vertex 
from the right can either come from mutually distinct 
or connected diagrams. The tree diagrams are a subset 
of the former group, and we argue that by letting Ai go 
to some bulk effective coupling Aoff all diagrams of the 
former group are included. The connected diagrams can 
be grouped by the number of times they are connected, 
and shown in Fig. ^(c) are a set of diagrams which are 
connected exactly once. Again we argue that by taking 
Ai — > Aefi the diagrams of Fig. ^(c) give the entire contri- 
bution of the set which are connected exactly once. The 
sum of these diagrams is evaluated in the appendix, and 
is found to contribute to A a term which is higher or- 
der in no than that given by the tree diagrams. It can 
be shown in general that the groups with more connec- 
tions will contribute correspondingly higher order terms, 
and therefore this classification scheme gives rise to an 
expansion for A. 

Of course, an almost identical mechanism will also gen- 
erate a term A'^^'"^^(0)^, and it is straightforward to show 
on the grounds of symmetry that A*^^'") = —A. Although 
this initial term is equally relevant from the renormaliza- 
tion group point of view, nevertheless it does not affect 
the late time behavior of the density. This is because it 
acts as a source for late time fluctuations only through 
the response function ( 2.20| ) , and this is strongly damped 
for t2 ^ ti. In contrast the response function of the ij: 
field is simply the diffusion propagator, which has no such 
damping. 

In summary of the discussion above, for 2 < d < 4 
and for large times one can replace the full theory with 
a simplified action 



S 



dt 



-Aeff0((/.2-^2)}-n^^(O)-^V;(O)^ 



-|-other initial terms 



(2.27) 




FIG. 5. Diagrammatic expansion for {</!>). Diagrams which 
contain initial terms other than and A are not shown, 
but are included in the sum. The only diagram in which the 
leftmost vertex is connected to 1/) fields is that of the single tp 
loop. 



terms within curly braces, is linear in and ip these fields 
can be integrated out to yield the equations of motion 



= V^^-Aeff^^-hAeffV'^ 

at 



(2.28) 



(2.29) 



These are equations for classical fields with fluctuations 
in the initial conditions. They are often taken to be 
the continuum limit of the master equation (2.1), but 
we stress that only for d > 2 and large times are these 
equations valid. In addition, since Acfi < A, it is never 
correct to say that (ab) ^ (a) (6) , but only that they are 
effectively proportional. 



III. DENSITY CALCULATION FOR EQUAL 
DIFFUSION CONSTANTS 

A. Effective action: 2 < d < 4 



Starting with the action (2.27) one can calculate ex- 
actly the leading time dependence of the density, as well 
as correlation functions. We begin with a comment about 
notation. For this section and the next, where we deal 
with only the effective field theory, averages over the ini- 
tial conditions will be denoted by angular brackets. The 
classical fields 0, tp themselves represent bulk averages, or 
equivalently, averages over reaction and diffusion, of the 
same fields as written in ( 2.18 ). Also, the effective cou- 
pling is abbreviated to be A — Acff . With this notation, 
then, the average of equation ( 2.28| ) over the translation- 
ally invariant initial conditions is 



dt 



(0) =-A(02)-(-A(7/;2^, 



(3.1) 



where A is given by (1.6). Since the bulk theory, the 



since V^icp) = 0. 

A diagrammatic expansion for (</>) is shown in Fig. 
Operating on both sides of this expansion with {dt — V 
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the inver se o f the Green's function propagator, gives 
equation (3J). At this point, knowing that A is relevant 
for d < 4, one might attempt to apply the renormaliza- 
tion group to try to find a nontrivial fixed point of order 
4 — d. However, no such fixed point exists. This is be- 
cause there are no corrections to a correlation function 
{ip{xi,ti)ilj{x2,t2)'ip{0)'^), so that A is not renormalized. 
It therefore flows, for d < 4, to infinity under renormal- 
ization. This would appear to make it very difficult to 
sum the diagrams in Fig. ^ explicitly. However, it turns 
out to be possible to solve ( ^T| ) exactly for late times. 
There is only one diagram contributing to the value of 
(tp^) in equation (3.1), which is the single i/' loop. Eval- 
uating this loop gives {ip'^) = A/(87rt)'^/^. It is impor- 
tant to note that this result holds even when all possible 
higher order initial terms are included. 

Next, consider a related problem in which in (3.1) 
is replaced by (0)^, which is equivalent to including only 
the diagrams in Fig. || which are disconnected to the 
right of the leftmost vertex. This partial sum satisfies a 
differential equation known as Ricatti's equation, which, 
though non-linear, can be solved.! Let / denote the func- 
tion which satisfies this equation, that is. 



-d/2 



(3.2) 



It will be shown below that this function / provides a up- 
per bound for the actual density, but first we will discuss 
the solution of this equation. It is integrable for certain 
values of d, specifically d = 4 and d = 4 ± 4/(2s + 1) 
where s is a non-negative integer. For general values of d 
a solution can be obtained by transforming the equation 
via the substitution / = u/{Xu), which gives 



A^A 

(8^ 



-d/2 



(3.3) 



a linear, second order equation whose solution can be ex- 
pressed in terms of confluent hypergeometric functions. 
Therefore the asymptotic behavior of / is rigorously ob- 
tained, and is in fact what one naive ly o btains by assum- 
ing / ~ and inserting it into (3.S): 



Ai/2(87rt)-''/4 d<4 
1 d = 4 

1 d > 4. 



Alt 



(3.4) 



[xty 



When d < 4 the asymptotic behavior comes from balanc- 
ing the two terms on the right hand side of ( ^.2| ) , whereas 
for d > 4 it comes from balancing the p and the / terms. 
For d = 4 all three terms contribute, and the amplitude 



^For an interesting mssentation of the properties and history 
of this equation, see |G1|. 



4" — — 



1 



(3.5) 



The case of d = 4 will be discussed in more detail in 
Notice that the asymptotic behavior of 



section [II B 



the solution / is independent of the initial conditions. 
In fact, the initial conditions must be specified at some 
to > 0, since the equation is singular at t = 0. A natural 
choice for this initial time is that given by the coarse- 
graining time scale of the effective initial conditions, that 
is to = (n^A)"^ 

Now we show that / provides an upper bound for the 
actual density (0). Our method is to derive an equation 
for X = / ~ (9^) s-nd show that asymptotically x > 0. 
Since (j> is a real field in the effective theory, then h{t) = 
{(fp) — (0)^ > 0. Equation (3.1) can be rewritten 



and then substituting {(j)) — f — X gives 



d 



-X = \h + \{x-2f)x. 



(3.6) 



(3.7) 



Assume that x(^o) = 0, that is we choose the initial con- 
dition for / such that /(to) = (0(io))- As mentioned 
above, the asymptotic value of / is independent of the 
choice of initial conditions. Since / is known to be pos- 
itive for all t > tg, then from equation (3.7) we know 
that X > whenever x < 0. Now suppose that there 
exists some ti > to for which xi^i) < 0. Since x(t) is a 
continuous function, it follows that there must be some 
intermediate time to < t < ti for which xit) < and 
x(t) < 0. This is in contradiction with equation (^), 
and therefore our assumption that there exists x(^i) < 
for ti > to is false. 

We can also find a lower bound for {(j>) by noting that 
0(x, t) > |'0(x, t)| at all points (x, t). This is equiva- 
lent to the statement that a(x, t), 6(x, t) are at all points 
non-negative, when starting from any initial condition in 
which a, b are everywhere non-negative. While this result 
is somewhat intuitive, it can be made more rigorous by 
considering the field equations (2.28), ( |2.29| ) expressed in 
terms oi a = {(j> + ip) / V2 and &==(</)- ijj) / V2: 



ot 



= V^b - V2Xab. 

at 



(3.8) 



Given that the fields a, b are initially everywhere non- 
negative, then for the fields to have a negative value at a 
later time ti there must be an intermediate time < to < 
ti for which both a(to) = and dta{to) < 0. However, 
in the case where a = at a single point in space, then 
a > locally around the point, implying that it is a local 
minimum and dta > 0. For a region of a = equation 
gives dta = in the region and dta > on the 



8 



boundary. Therefore the fields cannot pass through zero, 
and will remain non-negative. 

Since > it follows that {(/))> {\tp\). At late 
times Tp has a normal distribution, independent of the 
initial distribution, which follows f rom the fact that tp 
obeys the simple diffusion equation ( 2.29 ). Therefore the 
asymptotic value of (li/'l) can be computed directly. The 
asymptotic distribution of ip is given by 



P[ip{t)] cx exp <^ - 



mm 



from which it follows that 



(2A) 



1/2 



7^1/2 (87r)<i/4 



(3.9) 



(3.10) 



Given the upper bound {(j)) <f ^ 0{t-'^/^) it can be 
shown that {(p) ~ (|^|), that is, that the lower bound 
gives exactly the density. Since (g^) > (g)^ for any real 
g, then 



(^2)_2(0|V|) (3.11) 



Using again 



> 



{cp - ivi)^ < m - m = --(0) - o{t-'-'^'). (3.12) 

Therefore (0) = + 0{t-^^^-'^^^), which gives {(p) ~ 
(Itpl) for d < 4. This is actually a statement about segre- 
gation in the system, implying that to leading order the 
density of a -I- 6 is the same as |a — 6|, or equivalently, that 
the minority species in each region decays faster than the 
majority. For 2 < d < 4 then, we find 



(a) 



Ai/2 



^l/2(8^)d/4 



-d/4 



(3.13) 



as stated in section I, with A given by (1^). Substituting 



the leading order term in the expansion A — nQ + 0{nf/'^) 
then gives the result of Toussaint and Wilczek Q . In fact, 
our method is very similar to theirs, with two exceptions. 
First, they use a central limit argument to calculate A, 
whereas we can compute it directly from the full field 
theory. It is reassuring that the answers agree, to leading 
order in Uq. The other difference is that they calculate 
(IV'I), and then hypothesize the asymptotic segregation, 
s ayin g (0) ~ (IV'D- Starting from the effective theory 
(2.27) we have shown rigorously that these quantities are 
asymptotically the same. 



B. Effective Action: d > 4 

When d — 4 the upper and lower bounds for the density 
from section HI A still hold: (lipl) < (0) < /• However, 

(IV'I), since the 



it is no longer necessarily true that {4>) 
bound on the corrections, which is of order 0(i^^/^^'^/*), 



is the same order as the density. T he u pper bound (0) < 
f ^ A^/t is given by (|]|) and (U). Notice that for 
small A or small A that 1/A. Also, when A is 

large or A is large then A^/^/(87r). However, in 

the intermediate region there is a smooth crossover in the 
upper bound from the A dependent asymptote to the A 
dependent asymptote. 

The lower bound is given by (0) > — A\/t with 
— A^/^/Stt'^/^. For large A, then, the upper and 
lower bounds differ by a factor of ^/tt. The lower bound 
continues to decrease with A, and therefore is not very 
useful in the small A limit. However, since the parameter 
A is dimensionless in c? = 4 one can do a perturbative 
expansion for small A, which resu lts in a better lower 
bound. It follows from equation (3T) that the zeroth 



order term in this expansion is a constant, and is in fact 
equal to the small A limit of the upper bound, A^^. To 
the next order one has 



(3J4) 



and it is plausible to conjecture that the amplitude is 
monotonically increasing with A. 

The amplitude given by Bramson and Lebowitz Q , has 
the form 



A4 oc 



constant A < A„ 



Ai/2 



A > A, 



(3.15) 



that is, the amplitude is independent of A for small A. 
Their result seems to be at odds with our small A cal- 
culation. However, as discussed in section |l[ there are 
differences between our model and the one they study. 
Since the corrections to the small uq or A limit are non- 
universal, this is a possible explanation of the discrep- 
ancy. 

When d > 4 then it follows from the power counting 
of section 11 C that the (A/2)^^ initial term is irrele- 
vant. In this case the density is given asymptotically by 
(a) ~ {Xt)~^. The power law of the density decay is 
independent of the dimension of space. The amplitude 
A~^ will depend on the dimension and the microscopic 
details, but it is independent of initial terms, or equiva- 
lently initial conditions. 



C. Renormalization for d < 2 

When d < 2 one has to consider the full theory as given 
by the action (|2.18|), and the subsequent renormalization. 



Although (3.1) is still valid formally, since the noise in 
(^.23 ) averages to zero, we can no longer apply the upper 
and lower bounds of the previous section since, in the 
presence of the imaginary noise term, (p is no longer real. 

Much of the contents of this section is directly related 
to the one-species calculation of Ref. |l^ , in which more 
details can be found. The primitively divergent vertex 
functions were identified by power counting in section 
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X+t X 



- -ozk 



<:;o:: 



;< 



FIG. 6. The sum of diagrams which contribute to the prim- 
itively divergent vertex function Xi{k,t). Shown here is the 
case i = 1, with the ip propagators for the incoming external 
legs. 



[IB, and were found to be those with two lines com- 



ing in and two or fewer lines going out. These primitive 
divergences are used to define a renormalized coupling, 
following conventional RG met hods [2^ ] . It is found that 
all the vertices in the action (2.18) renormalize identi- 
cally, with the primitive divergences given by the bubble 
sums shown in Fig. ^. 

In this sum all diagrams of a given number of loops 
come in with the same sign, since replacing a (j) loop with 
a ip loop, for example, introduces always two negative 
signs (see Fig. |l]). At the order of n loops there are 2" 
diagrams, so these form a geometric sum, with the ratio 
given by 2 times the value of a single loop. Denoting 
this sum by Xi(k,t) where i = 1,2 labels the number 
of outgoing lines, then the Laplace transform, Xi(k, s) — 
/(J dt e~^*Ai(fc, t) is given by 



Ai(fc,s) 



1 + [A/iSirY/^] \2 r(e/2)(s + fcV2) ' 



(3.16) 



where e = 2 — d. The renormalized coupling is defined 
in terms of an arbitrary normalization scale k, which has 
dimensions of wave number: = K~'^\2{k = 0, s = k^). 
Then the (3 function is 



OK 



(3.17) 



which gives a fixed point .g|j. = 0(e). 

Let the density n{t) = {a{t)) = {b{t)). Since the 
density is independent of the normalization scale, then 
dn/dn — 0, which leads to the Callan-Symanzik equa- 
tion 



where in the asymptotic limit of large K^t the running 
coupling has the limit gji — > g^. However, the run- 
ning initial couplings go as = n^{K'^t)'^^'^ and A — 
A(K^i)^/'^, that is, they flow to a strong coupling limit. 

The solution ( |3.19 ) is used to calculate the asymptotic 
density in the following way: the density is calculated as 
an expansion in gji and n^, a nd this expansion is put into 
the right-hand side of (3.19). Then, in the limit of large 



= 0. 



n{t,gR,n^,l\) 
(3.18) 



The solution is found by the method of characteristics to 
be 

n{t,gR,n^,/^) = {nH)-^/^n{K-^,gR,fi^,A), (3.19) 



K t, the coupling expansion will yield an e expansion, 
but only if the behavior at large ri^ and A is controlled. 
This may be done if the diagrams may be grouped into 
sums over all powers of and A, which, when summed, 
yield a well-defined limit. In the one-species case this 
grouping was relatively simple. The series may be put 
in the form of a sum of terms gfl?^0/n(ffi^"■0), where n 
counts the number of loops in a given diagram. The 
term n = corresponds to the sum of tree diagrams, 
given by the solution of the simple rate equation, so that 
foigRntj,) oc {gjih^)~^ as oo. By explicit calcula- 

tion, it is then possible to show that the /„ for n > 
behave in a similar manner. Since gn g^ = 0(e), 
this lead to the result that n{t) ~ Ajf^l"^ where the am- 
plitude A is in principle calculable to any order in e. 
In the present case, the series may similarly be orga- 
nized as a sum of terms of the form 5^?i0/n((?i?«-0, <?|fA)i 
where now the n = term is given by the sum of dia- 
grams in Fig. H This is given by the solution of ([3.l|), 
which, by the analysis of the previous section, implies 
that /o(5i?n0,g|jA) cx {gj^Ky/'^ / {gnn^) in this limit. 
However, unlike the single species case, n does not sim- 
ply count the loops, since already at n = there are 
arbitrarily many ■0 loops. In addition, while it is pos- 
sible to identify those diagrams appearing at ti = 1 for 
example, it is difficult to see how to ex pre ss their sum in 
terms of a suitable generalization of (3^), and thereby 
evaluate it. Assuming that their asymptotic behavior is 
independent of no and thus of n^, there are three conceiv- 
able ways in which these higher order terms could affect 
the result. They either diverge less slowly than A^/^ as 
A — > CX3, in which case the n = result gives the leading 
behavior, which would then yield the same result as for 
d > 2; or they all behave like A^/^, in which case the den- 
sity behaves as i"''/^ but with an amplitude which has a 
nontrivial expansion in powers of e; or they diverge more 
strongly, in which case the density no longer behaves as 
t~'^/^ for d < 2. Since this last possibility is in conflict 
with numerical experiments and rigorous results (albeit 
for slightly different models), it is unlikely to be correct. 

When d = 2 the running coupling goes to zero as 
{lu.t)~^ for t ^ oo, rather than to an order e fixed point. 
Therefore the leading order terms for an e expansion of 
the amplitude become the exact asymptotic amplitude, 
with correction terms which go as (Ini)^^. Therefore, if 
our conject ure i s correct, then density should be given 
exactly by (3.13) in the large t limit. 
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IV. DENSITY CALCULATION FOR UNEQUAL 
DIFFUSION CONSTANTS 



When the two species of particles no longer have equal 
diffusion constants, then the vertices which depend on 5 
must be included in the full theory. Then for d > 2 an 
effective theory can be developed, just as before, with the 
resulting action 



H 



5 = d^x 



dt 



[Kdt 



A 



^ ^ t t , 



(4.1) 



The effective theory describes classical fields which evolve 
via the deterministic equations of motion 



d 



(4.2) 



(a) 



(b) 



(4.3) 



which follows from integrating out the 0, ■)/) degrees of 
freedom in the bulk component of ( 4. 1 ) . From these equa- 
tions the density can be calculated exac tly b y using the 
same methods as before. First, equation ( 4.2) is averaged 
over the initial conditions to yield equation (3.1), just as 
in the (5 = case. The solution to Ricatti's equation 
again provides an upper bound (0) < f ^ \/ {ijj{tY), al- 
though the value of {ip{t)^) is changed. It will be shown 



that (V'^) oc r 



-d/2 



SO the upper bound decays with the 



same exponent as before. Since the fields are real and 
4> > IV'I, it then follows that {(j)) ~ for d < 4, 

as shown in (3.12). Furthermore, it will be shown that 



asymptotically ipit) has a normal distribution, so the den- 
sity is given exactly by (a) = {(j)) /V^ ^JJ^P^YpK. There- 
fore the only change in the asymptotic density from the 
(5 = case is due to the change in the value of (-0(i)^). 



A. Calculation of i,i){tf) 

The initial terms in the effective theory are in gen- 
eral changed by the presence of 8 in the full theory, and 
therefore must be computed again. One can show that, 
as before, no A^^'^^^/i initial term is generated. For any 
diagram which ends with a single "0 line, the last vertex 
(first from the left) must be a vertex. However, this 
external line has fc = 0, and so all of these diagrams have 
no contribution. To leading order A = no is unchanged, 
as can be seen from the diagrams in Fig. ^: the leading 
order contribution to A comes from diagrams composed 
of no loops, and so all lines carry wave number = and 
are unaffected by the (5fc2 vertex. The correction terms 



FIG. 7. The response functions for the case 5 7^ 0, and the 
coupled equations they satisfy. 

to the small no limit of A will likely be of the same order 
as before, Oiv^J'^X^I^)^ but with a different amplitude. 
This amplitude could be calculated, although it would re- 
quire a generalization of the response functions discussed 
below. It will be shown the asymptotic value of (V)^) de- 
pends only on A, and so the other surface terms can be 
neglected. 

There are new response functions generated in the bulk 
theory. With (5 = there was just a bare -0 propaga- 
tor and a response function. In this theory there are 
instead four response functions, which connect 0,';/' to 
(^, in each possible way, as shown in Fig. ^ Each of 
these response functions, represented by double lines, is 
an infinite sum over all possible numbers of (5^^ vertices 
inserted. 

These response functions can be found exactly via cou- 
pled integral equations, also shown in Fig. |^. For our pur- 
poses, since just the leading term for small uq is being cal- 
culated, we need only the form of the response functions 
when the earlier time argument is set to zero. To calcu- 
late the higher order terms in the expansion A = no -f . . . 
one needs to derive these response functions with ti ^ 0. 
Setting ti = t\ — Q in the equations represented by 
diagrams (a) gives 



dt'e-'''''^*-'">F{k,t') (4.4) 



F{k,t)^5e / dt'e'^"^'-'"^ 



G{k,t') 



(4.5) 
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or, in terms of f,g defined by G{k,t) = e ^^^g{k,t) and 
F(fc,t) = e-'='*/(fc,i) 



g{k,t) ^l + 5k^ / dt'f{k,t') 
Jo 



(4.6) 



Differentiating equations ( |4.6| ) and ( [4.7[ ) with respect to 
i gives 



fik,t) = ^gik.t) 



(4.8) 



— [(1 + n^Xtffik, t)] = <5fc2(l + n^Xtfgik, t) (4.9) 

Substituting for / into the lower equation and manipu- 
lating the expression gives a remarkably simple equation 
for g 



o2 

— [{l + n^\t)g]^5^k%\ + n^\t)g] 

which has the general solution 
1 



(4.10) 



g{k,t) = 



1 + 710 At 



cosh(<5fc2t) + Bsinh((5fc2t)] 



(4.11) 



From the integral equation (4^) one finds the conditions 
(^(/e, 0) = 1, which implies A — and g{Q,t) — 1, which 
then implies B = np A/ ((5 fc^). Therefore the explicit form 
of G{k,t), and from (|4.8|) F{k,t), is calculated: 



G(/fc,t) = 



F(fc,i) 



1 + n^At 



(1 + n^xty^ 



coshiSkh) + ^ smh{Skh) 



Sk^ 



(4.12) 



1 + n^At - ^Itx I sinh((5fc^i) 



■ cosh{SkH) 



(4.13) 



The other response functions, H{k,t) and J{k,t), de- 
fined in diagram Fig. 0, can be found via similar meth- 
ods. The coupled integral equations shown in Fig. |^(b), 
written in terms oi h — *H and j = * J, are 



h{k,t) 



{1 + n^Xty 



Sk^ f dt'f- 

Jq V 1 + '^'t'Xt J 



(4.14) 



J A 



+ irrelevant terms 



FIG. 8. The generalization of the simple ip loop of Fig. ^ 
to the case of 5 7^ 0. 



j{k,t) = Sk^ / h{k,t') 



(4.15) 



Differentiating both equations with respect to t and sub- 
stituting to eliminate h gives the equation 



[{1 + n^Xt)j] = d'k%l + n^Xt)j] 



(4.16) 



which has the general solution 
1 



1 + n^Xt 



[Acosh{Skh) + Bsmh{Skh)] 



(4.17) 



The condition that j(fc,0) = impli es A = 0. The gen- 
eral solution of h can be found from ( 4.17 ), and then the 
condition that h{k,0) = 1 implies B = 1. Therefore H 
and J are given by 



H{k,t) 



{l + n^Xtf 



{1 + n^Xt) cosh{5k'^t) 
nj,X 



sinh((5fc^i) 



J(fc,t) 



(l + n^Xt) 



smh{Skh). 



(4.18) 



(4.19) 



In section III A the value of UjP) was calculated from 
the simple loop shown in Fig. The generalization of 
this calculation is given by the diagrams shown in Fig. 
which are composed of the G{k,t) and J{k,t) response 
functions. The surface couplings A^^'^) ^ —/^i^fi) 
yond the leading small no terms, and so the couplings are 
labeled A and A' respectively. It should be noted that 
unlike the ^ = case, these are not the only diagrams 
which contribute to (V'^)- Examples of other diagrams, 
and arguments for why they are irrelevant, will be given 
below. First, we compute those of Fig. H, which give 



_c£k_ 
(2^ 



[AG(fc,t)2- A'J(fc,t)2]. (4.20) 



Substituting ( 4.12 ) and ( 4.1£ ) into the equation above, 
and rewriting the integral in terms of the variable u = k'^t 
gives 
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1.5 



1/2 



d=l 
d=2 
d=3 
d=4 



(a) 



(b) 



A.'k^\(/(|) 2 




FIG. 10. (a) An example of one of the diagrams besides 
those of Fig. ^ which contribute to {ij)'^), and (b) the effective 
bulk vertices that all such diagrams contain. 



FIG. 9. A plot of ^/Q = {a)s/{a)o for integer values of d. 



(4^)'i/2r(d/2)(l + n^At)2 7o 



Acosh"'((5u) - A' sinh2((5u) 



(5w 



sinh((5w)^ 



sinh(2(5w) 
(4.21) 



(5m 



Each term in the square brackets gives a convergent in- 
tegral for c? > 0. Therefore we can take the large t limit 
before integrating, and only calculate the leading term 
in t, which is found to be that on the far right in the 
brackets. Consequently, the value of A' is unimportant. 
Evaluating this integral gives 



A 



:Qid,d) r^/' 



where 



Q{d,6) = 



4 [(1 + 5)2-'^/2 + (1 - 5)2-<i/2 _ 2] 

,52(d-2)(d-4) 



(4.22) 



(4.23) 



From ( |4.22| ) it follows that Q ^ (V'^)5/(V'^)o, in the small 
Uq limit. This function Q is non-singular at S — 0, and 
satisfies Q{d,0) = 1. While Q appears to be divergent 
at d = 2,4, it is actually finite everywhere except d > 4 
and (5 = ±1. It is likely that the limits of t ^ oo and 
(5 — !■ ±1 do not necessarily commute, and that a separate 
treatment for the case of an immobile species, at least 
in this singular case, would be required. For d < 4 this 
function has finite values as S ±1, but the slope at 
6 = ±1 is infinite for d>2. 

While the calculation of Q{d,6) is only strictly valid 
for 2 < d < 4, it is nonetheless interesting to consider its 
limits for the integer dimensions from d = 1 to d = 4, 
motivated by section IIIC on d < 2, in which it was con- 



jectured that the "classical" amplitude is als o the leading 
term in an e expansion for d = 2 — e. From (4.23) 



Q{d,6) 



^[(l+5)3/2 + (l-<5)3/2_2] 

(1-^) \n(l - 5) + {1 + S) \n{l + 5) 



j^[2 - VT+s - VT^] 



d = 1 

d = 2 

d = 3 

d = 4 
(4.24) 



Since the density goes as ("0^)1 the function Q{d, 6) 
is plotted in Fig. ^ for integer values of d. The density 
amplitude increases monotonically with \d\, but is not 
changed remarkably for modest values of S. 

There are other diagrams which give contributions to 
(-0^), unlike the ^ = case. Some of these are shown in 
Fig. |l^(a). All of these diagrams have the similar fea- 
ture that they contain one of the two sub-diagrams in 
Fig. |l^(b). These sub-diagrams give rise to effective ver- 
tices of the form AVV202 and A'V'V2(/)2 in the bulk the- 
ory. However, such vertices are irrelevant, which follows 
from power counting, and so the diagrams which arise 
from them must be sub-leading in time. Therefore we 
con clude that asym ptotically the value of {ip'^) is given 
by (imi) and ( ^!23| ). 



B. Demonstration that tp{t) has a Normal 
Distribution 

In order for the calculation of {ip"^) to give the ampli- 
tude of the density it is necessary that ipit) have a normal 
distribution. When S — this follows directly from the 
simple diffusion equation satisfied by -0, or equivalently, 
from central limit arguments. However, ^ evolves via 
equation ( [1.3| ) for (5 7^ 0, and so it needs to be shown 
that it still flows to a normal distribution. What we will 
show is that the random variable i'*/^?/' flows to a static 
normal distribution, the width of which was calculated 
above. 

Consider where n is even. There is one diagram 

in which n response functions G{k,t) are connected in 
pairs to n/2 initial terms (A/2)'02. This diagram con- 
tains n/2 loops, and is therefore of order It 



13 



was shown above replacing any of the G{k,t) loops with 
J{k,t) response functions connected to (A'/2)0^ gives a 
lower order contribution. Similarly, any other diagrams, 
which would originate from considering higher order sur- 
face terms, will involve more than n/2 loops, and will 
therefore decay faster in time. For n odd one finds that 
there are no diagrams for ("0") which decay as slowly as 
That is, for n odd, limt_oo((i''/'*V')") = 0. Since 
the distribution of the variable f^'^'^i/j has only even mo- 
ments as i — > oo, and these moments are just multiples 
of {{f^/^ijjY), generated by all possible pair contractions, 
then the distribution is normal. 



V. CORRELATION FUNCTIONS FOR 2 < d < 4 

When d > 2, one can use the classical action to cal- 
culate the correlation functions. Consider the distribu- 
tion of the ra ndom variable i''/*(/)(x, t) with 2 < d < 4. 
we know that {f^^-^cj) - 



From section III A we know that - t'^/'^lipl) 

as i ^ oo. Furthermore, from equation ( 3.12 ) it follows 
that, as t ^ oo, {{t'^/'^(f> - f^/^li;])^) 0. This suggests 
that the distributions Plf^^^cj)] P'lf^^^li/j]] as i -> oo. 
The latter distribution is known exactly, as f^^"^^ is at 
late times given by a static normal distribution. 

It is not correct to say that asymptotically (p and 
are everywhere equal, since this would imply that there 
are no regions in which the densities a and b are both 
non-zero. However, the reaction regions, those in which 
both densities are non-zero, become negligibly small for 
large t, and the corrections to setting </> equal to jV'l in cal- 
culating correlation functions will be subleading in time. 
Stated another way, the leading term in both (0i02) and 
(|V'i||V'2|) is of order t^'^^^. To this order the two random 
variables (j) and |V'| have identical distributions. This is in 
contrast to a quantity such as cfP' — ijP', which is measuring 
a subleading term relative to 

We can use the property that f^/^cj) is given by the 
absolute value of a gaussian random field to calculate 
correlation functions. This is similar to what is done is 
the dynamics of phase ordering, where the order param- 
eter field can be mapped to an auxiliary field which is 
assumed to be a gaussian random field. This analogy 
will be discussed further below. 

Since (j) and |^/^| are given by the same distribution, 
we conclude {4'i4'2) ~ (|V'i||'02|), where the labels indi- 
cate the positions Xi and X2 at time t. The correlation 
function (|V'i||V'2|) can be calculated exactly by using the 
fact that, asymptotically, ip{t) has a normal distribution. 
The joint probability distribution P[?/'i,-02] is then also 
normal, so 



= ^ ^ exp{-a0i^-aV^|-/3ViV'2}, 



27r 



(5.1) 



where we have used translational invariance to set (-01) = 
(?/^|). The constants a and (3 are determined by the 



<\i/(k)\|f(-k)> = ^ A 

FIG. 11. The diagram for (V'(k)V'(-k)), when 5 = 0. 



values of (0^) and (V'i'02), which are evaluated from 
the diagrams. The latter we have only calculated for 
i5 = 0, or equal diffusion constants, so we consider 
that case first. For notational convenience we define 



C{t) = A/(87rt)''/2. The diagram shown in 



Fig. 0(a) is used to calculate the correlation function 
('0(k)'0(— k)), from which one finds 



(2^' 



ik-(xi— X2) 



(0(k)VX-k)). (5.2) 



When (5 = then (V'(k)V'(-k)) = Ag-^'^"*, and 

{M2) = Cit) exp(-rV80 = C{t)f{r^/t) (5.3) 



where r = |xi — X2I. In terms of (5^) we find for a, (3 
1 ^ / 



2C{\-p) 



c{i-py 



(5.4) 



With these values substituted into (|5.lD, one can cal- 
culate 



(0102) - (|0'l||^2|) 



2C 

TT 



\/l — P + f arctan 



d^Ji I d4,2\MHP[^i.H 
f 



(5.5) 



This correlation function can be used to find the correla- 
tion functions (aia2) and (0162)- Specifically 



1 



(0102) 



102 + '010'2), 



(5.6) 



which gives for the connected part (0102)0 = (01O2) - 

(«>^ 

(aia2)c = 

, -/- 1 + Vl - + /arctan I JL^ ] 



For large r, f = exp(— r^/8i) is small, giving 



(aia2)c 



2{8Trt)<if^ 

Similarly, (0162) = (0i02 ^ 0i0'2)/2, so that 



(5.7) 



(5.8) 
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<aa> 



0.4 




-0.4 



FIG. 12. The correlation functions 
{a{r,t)b{0,t))c plotted as functions of r 
is given in units of A(87rf)~'^''^. 



a{r,t)a{0,i))c and 
t. The vertical axis 



{aib2)c 
A 



/ 



7r(87rt)''/2 

which for large r goes as 
(ai&2 



2 / - 1 + VW^ + / arctan ( 



A 



4(87rt)''/2 



(5.9) 



(5.10) 



A plot of these connected correlation functions is shown 
in Fig. The signs (0102) > and (0162) < can be 
understood for short distances to be a consequence of the 
segregation. Given an A particle at a particular point, 
there is an increased probability that a nearby particle is 
also an A, and a decreased probability that it is a B. 

F or th e case 5 7^ one has {iP) — C{t)Q{d, S), as given 
by ( 4.22 ). The generalization of ('!/'(k)?/;(— k)), shown in 
Fig. 11, behaves for small k the same as when 6 = 0. 
Therefore, for large r one still has {1P11P2) = Cf. When 
this is put in the expressions for (0102) and (0162} one 



finds that the large r behavior is given by (U) and ( |5T0| ) 
is unaffected by 5 ^ 0. 

While these correlation functions and other quantities 
can be calculated, they ultimately rely on the stronger 
statement that '0 is a gaussian random field, and that 
~ \^p\. The topology of the domains is determined by 
the random field, with the boundaries between a regions 
and b regions given by the zeroes of ip. This topology is 
completely equivalent to an analogous situation in phase 
ordering. It has been suggested that in the phase order- 
ing of a scalar order parameter an invertible, non-linear 
mapping from the order parameter field to an auxiliary 
field results in the latter being a gaussian random field 
p3[ . Usually this mapping is chosen to be the solution of 
a single kink, for example the hyperbolic tangent profile. 
While this method is no longer believed to be quantita- 
tively correct EJ] , it does provide a useful picture of the 



structure of the domains. Again, the zeroes of this gaus- 
sian random field determine the boundaries between the 
equilibrated phase. 

The difference between these systems lies in how cor- 
relation functions are calculated from the random field. 
In the reaction-diffusion case one is interested in the cor- 
relation functions of the field itself, and of the absolute 
value of the field. Neither of these quantities exhibit re- 
markable behavior. In the phase ordering one argues that 
at late times the mapping between the order parameter 
field and the gaussian field goes to a step function, and 
therefore order parameter correlations are given by the 
correlations of the sign of the random field. These sharp 
boundaries give rise to more interesting features, such 
as non-analytic terms in the small r limit of the correla- 
tion function, or correspondingly power law tails for large 
wave number in the Fourier transform. 



VI. REACTION ZONES 



It was shown in section [II A that for d < 4 the par- 
ticles segregate asymptotically into regions of purely A 
or B particles. As a result of this segregation there exist 
interfaces between the two species, and all reactions oc- 
cur in the interfacial regions. These reaction zones have 
interesting scaling properties. For example, the width 
of the interface goes as with the exponent a < 1/2. 
Also, the nearest neighbor distance distribution of the 
particles in the reaction zone is found to have a charac- 
teristic length £rz that goes as a power of t, with an ex- 
ponent which differs from that of the bulk system, where 
(a)"^/'' ^ t^/^. To derive these properties we begin with 
a related steady-state problem. 

Consider a system with a source of A particles located 
at the boundary x = —L which maintains a fixed current 
Jx, and a similar source — Jx of B particles positioned at 
X = L. These opposing currents will establish a steady- 
state profile, in which the average densities will be func- 
tions of the transverse coordinate x. For a given current 
J one can choose L to be large enough that the reactions 
are localized to an interfacial region of width w L. 
In this case, it is found that the densities in the reac- 



tion zone, where 



< 



have universal scaling forms. 



Also of interest is the reaction rate R{x) = Xo{a{x)b{x)), 
which exhibits scaling, and which is used to define the 
width of the reaction zone. 



The power counting of section 11 B showed that the 
four-point vertices were irrelevant for d > 2. Therefore 
R{x) ~ Aoff(a)(6) in the asymptotic limit — which will be 
shown to be the small J limit — and the problem reduces 
to the differential equations of the effective theory: 

^ ' - = V2(a> -Aeff(a)(6) (6.1) 



dt 
d_ 

di 



(&) = V2(6) -Aeff(a)(5). 



(6.2) 



From these equations it has been shown that 



15 



(6.3) 



implicitly by Galfi and Racz [g^ , and lat er explicitly by 
Ben- Nairn and Redner From ( |6.3D one identifies 

the width w ^ J^^^"^, and the characteristic length of 
the particle distribution within the reaction zone £rz = 
{a{x = 0))"^/'^ J"^/^'^. The latter quantity is derived 
in Ref. 

For d < 2 one does not have simply differential equa- 
tions, and the full field theory must be taken into account. 
We begin by observing that the current Ja is given by 
a*dxa in the notation of section [IB , and similarly for 
Jb- From dimensional analysis [J] — fc''+^. 

We proceed with the renor malization of the theory, as 
was sketched in section III C . A normalization scale k is 



introduced, and used to define the renormalized coupling 
gu- Since physical quantities, such as the width, cannot 
depend on k, then 



K — W 

an 



d_ 



- /3(gi?,) 



d 



wiJ,gR,S,K)^0. (6.4) 



Note that, since there ar e no diagrams which can dress 
the two-point vertices in (2.18), 6 does not get renormal- 
ized, and therefore does not appear in equation above. 
From dimensional analysis one has 



d_ 



+ (d+l)J^ + l 



w{J,9R,S,K.)^0. (6.5) 



Combining these equations gives the Callan-Symanzik 
equation 



w = 0. 



with the solution 

wiJ,gR,6) = >iJ-'/^'+'^wiK'+\gH,S). 



(6.6) 



(6.7) 



In the small J limit then gji — s- 5]^, and the right-hand 
side is given by 



(6.8) 



Following the same procedure for any dimensionful quan- 
tity results in the scaling behavior being given by dimen- 
sional analysis. That is, £rz ^ w, 



and 



(a), (6) ^ J'^/('*+i)i^,^b(a; ji/('i+i)), (6.9) 



Rix) ^ j('^+2)/('^+i)G(a; (6.10) 



Note that these results imply that 
R ^ j(2-d)/(<^+i)(a)(&). This can be shown exphcitly by 
calculating R oc {ab) = A(J)(a)(6), where A(J) is given 
by the bubble sum shown in Fig. 0, with s = J^/ ("^+1) ^ k = 
0. In the smah J limit then A( J) - ji^-d)/{d+i) _ 



Since the current J may be thought of as being due to 
localized sources of A and B particles at x = ±L respec- 
tively, the coupling constant power counting arguments 
are formally the same as those of Ref. |l^ (see Sec. Ill C ) , 
in which the sources are localized at t = 0. |2^] Thus the 
various scaling functions above may, in principle, be cal- 
culated as an expansion in g^ — 0(e), in which the lead- 
ing term is given by the solution of the rate equations 
( |6.l[ ). The next order corrections to the reaction profile 
have been computed in Ref. |^ , where it was shown that 
the fluctuation corrections lead to a universal power law 
tail in this function. 

For d — 2 one has ga ~ 1/| In J| for small J, and the 
leading order result is therefore found by substit utin g 
this behavior into the solution of the rate equation (|6.l|). 
This leads to the resultsi 



w cx 



In J| 



1/3 



R{x) - A(J)(a)(fe) 

^ j4/3|lnjri/3G(^ji/3|lnJ|-i/3)^ 

4zOC J-i/3|lnJ|-i/6. 



(6.11) 

(6.12) 
(6.13) 



As was discussed in Ref. , the corresponding results 
for the time-dependent cases of segregated initial condi- 
tions or of randomly homogeneous initial conditions is 
given by substituting J (x or J c>c t'^(d+2)/4 j^-^j^j^ 

d = 2), respectively, in the above formulas. These results, 
for the case of segregated initial conditions, have been 
obtained recently via heuristic arguments by Krapivsky 
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APPENDIX: CALCULATION OF THE LEADING 
CORRECTION TERM FOR A 

In order to calculate the first order correction term 
to the expansion A = Aq -|- . . . we must first comment 
on the bulk diagrams which generate Aeff- The effective 
coupling can be calculated as an expansion in the bare 
couplings, via the diagrams shown in Fig. 13, The loop 



integrals in this expansion require the cutoff A, and one 
finds 



* These logarithms were incorrectly omitted from Ref. p^ . 
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k=0 



-^2 ^2 



FIG. 13. The expansion for the eflective couphng constant. 
The wave number integrals are regulated by a cutoff A. 



4A 



d-2 



(87r)''/2(d-2) 



OiXl). 



(Al) 



If the response functions in the loop of Fig. ^(c) were 
instead just propagators, then this set of diagrams would 
be included into those of F ig, ^(b) when the substitution 
A2 — > Aoff is made via (Al). Therefore, the terms which 
are new and constitute a correction to Aq are those in 
Fig. ^(c) with the propagator loop subtracted out. We 
define the large t limit of these diagrams to be Ai, that 
is 



Al 



dt2 I at I 
Jo 



-2k^t 



(1 + n0Aiti)2 



(1 + n^Aita)^ 



(87r)'i/2(d-2)(l + n^Ait2)2 



(A2) 



Performing the wave number integral with the A cut- 
off imposed in the same manner as in ( |AlD gives 



Al = 



dt2 



{SttY/^ Jo (l + n^Ait2)4 

(l + n^Aiii)2 Ad-^ 



dti 



{t2 



d-2 



(1 + n0Aii2)^ 



(A3) 



The ti integral can be evaluated as a Laplace convolu- 
tion integral, and the cutoff dependent terms cancel. The 
remaining t2 integral is 



Al = 



dt2t\""^ 



(87r)'i/2(d-2)7o (l + n0Ait2)4 

4/10^2 , 



1 



A-d {A-d){Q~d) 
This integral can be done exactly, giving 



(A4) 
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Al 



^2 



(d + 2)(rf + 4) 

1 sin(7r(i/2) 



(A5) 



In terms of the initial density n p an d the effective cou- 
pling then one finds the result (|l.6|) for A. Evaluating 
the diagrams such as those in Fig. ^(c), but containing 
more loops will then give the higher order terms in this 
small no expansion of A. 
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